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Abstract: In this paper wc describe how three qubit entanglement can be an- 
alyzed with local measurements. For this purpose we decompose entanglement 
witnesses into operators which can be measured locally. Our decompositions are 
optimized in the number of measurement settings needed for the measurement 
of one witness. Our method allows to detect true threepartite entanglement and 
especially GHZ-states with only four measurement settings. 

1 Introduction 

Entanglement is one of the most puzzling features of quantum theory and of 
great importance for quantum information theory. It is the resource that makes 
various quantum protocols possible that perform certain tasks better than it 
would be possible with purely classical methods p^. The investigation and 
characterization of entanglement with experimental and theoretical methods 
itself is therefore a task of great importance. 

From an experimental point of view it is important to find effective tech- 
niques for the production and detection of entanglement. For the detection of 
entanglement several strategies are known: One can determine the whole density 
matrix |5J and then try to apply one of the necessary or sufficient entanglement 
criteria, e.g. the PPT criterion One can also look for a violation of Bell 

inequalities 0. Furthermore, there have been several proposals of detecting 
entanglement without estimating the whole density matrix 

However, all these nice ideas have also some disadvantages. Estimating 
the whole density matrix requires a lot of measurements, and in fact one does 
not need to know the whole matrix in order to check whether it is entangled 
or not. On the other hand, looking for a violation of Bell inequalities may 
not suffice for a making a decision, since there are entangled states which do 
not violate any known Bell inequality |SJ E3] ■ It has even been conjectured 
that entangled states with a positive partial transpose do not violate any Bell 
inequality at all . Finally, the recent proposals of detection require collective 
measurements on several qubits which makes them hard to implement with the 
present techniques. 
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Furthermore, the schemes mentioned above are in some sense too general 
for many experimental situations. They assume that no a priori information 
about the state is given. However, in a typical experimental situation one often 
tries to prepare a certain pure state. Due to imperfections of the experimental 
apparatus the output state will be a mixture of the desired state and some noise 
instead. In this case it is desirable to know whether the produced state is still 
entangled or not. 

In |12l I13j we proposed a general scheme for entanglement detection for the 
case that some knowledge about the state is given. In this scheme we only want 
to use local projective measurements, since these measurements can easily be 
implemented in a lab. In addition, we would like to decrease the number of 
measurements needed to the minimum, of course. 

The scheme relies on the well known concept of witness operators An 
hcrmitean operator W is called a witness operator (or entanglement witness) 
detecting the entangled state g e if Tr(Wg e ) < while Tr(Wg s ) > holds for 
all separable states g s . Thus if a measurement yields Tr(W g) < then the state 
g is entangled with certainty. As a consequence of the Hahn-Banach theorem 
it follows that for every entangled state g e there exists such an entanglement 
witness and for many states it is known how to construct such witnesses |15j . 
After having constructed the witness, we decompose it into a sum of local mea- 
surements, then the expectation value can be measured with simple methods. 
This decomposition has to be optimized in a certain way since we want to use 
the smallest number of measurements possible. 

Our paper is organized as follows: The first section deals with the decom- 
position into local measurements because this is the core of our approach. We 
define there what we understand by an optimized decomposition. We would like 
to remark here that finding optimized decompositions is in general a hard task, 
much harder than constructing entanglement witnesses. In the second section 
we illustrate our approach with an example of a two qubit system. We construct 
the witness and determine its optimal decomposition. Finally in the third and 
main part we apply this idea to three qubit systems. We explain how GHZ-type 
and W-type entanglement can be detected with local measurements. We also 
determine the minimal number of measurements needed for this. 

2 Local decompositions 

Assume that we have an hermitcan operator H acting on a tensor product 
H = Ha®T~(-b <£> ■■■ ®'Hz °f two or more (but finite) dimensional Hilbert spaces. 
In order to slenderize the notation we look here at the case that we have a 
bipartite N x TV-system: H = Ha ®Hb with dim(W^) = dim(W B ) = N. But 
all definitions in this section can be extended to more parties in an obvious 
manner. In order to measure the expectation value of this operator locally, we 
have to decompose it in a sum of tensor products of operators acting on only 
one system, or, cquivalently, we have to decompose it into a sum of projectors 
onto product vectors: 

m n 

H = Ai ® B t = *\ e i)( e i\ ® l/<X/i|- (!) 

i=l i=l 
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By measuring the expectation value of the projectors | e*) (e* | ® | fi) (fi | and adding 
the results with the weights Cj this decomposition JIJ can be measured locally. 
There are, of course, many possibilities of finding a decomposition like So 
we have to optimize the decomposition in a certain sense. But there are even 
several possibilities of defining an optimized decomposition. 

On the first sight one might try to minimize the number of product vectors 
corresponding to minimizing n in Q. This optimization procedure is already 
known from the literature, it was considered in |16j . There it was shown that 
for a general operator acting on two qubits one needs five product vectors and 
also a constructive way of computing these product vectors was given. 

However, since we want to construct an experimentally accessible scheme 
for entanglement detection it is natural to look for a decomposition where Al- 
ice and Bob have to perform the smallest number of measurements possible. 
By "measurements" we understand here von Neumann (or projective) measure- 
ments, since they can be easily implemented. Such a measurement for Alice 
corresponds to a choice of an orthonormal basis of Ha , and Bob has to choose 
an orthonormal basis TIb, too. So any operator of the form 

N 

M= c M \e k ){e k \®\f l ){f l \ (2) 

k,l=l 

with (e s |et) = (f s \ft) = fist can be measured with only one collective setting 
of measurement devices of Alice and Bob. Alice and Bob can distinguish the 
states \ekfi), measure the probabilities of these states and add their results with 
the weights cm to measure M. We call an operator which can be measured with 
one measurement setting (like M in Eq. (j2J) a local von Neumann measurement 
(LvNM). 

Having understood what can be realized with one measurement setting, we 
can state another optimization strategy. We want to find a decomposition of 
the form 

k N 

F = EE 4l4)(4l®l//)(/ii (3) 

i=l k,l=l 

with (e\\e\) = (fl\fl) = fist and a minimal k. This k is the number of collective 
measurement settings Alice and Bob have to use in order to measure H. This 
optimization strategy is the aim we are considering in this paper when we talk 
about "optimized" decompositions. 

The reader should note that minimizing m in Q is not the same as our 
optimization strategy. However, for systems of dimension N greater than 2 it 
might be more suitable to decompose the witness as W = Y^T -A* <g> Bi . As 
shown in |17j , the expectation values of operators Ai or Bi can be measured by 
a POVM with a single qubit as ancilla instead of counting clicks for all possible 
N outcomes of the operator. Also minimizing the number of product vectors 
{i.e. minimizing n in is not the same. This will become clear in a few 
seconds, when we study two qubits. 

3 Two qubits 

We illustrate the method by considering an experiment that aims at producing a 
certain 2-qubit state |\E') = a|01) + 6 1 10) written in the Schmidt decomposition, 
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i.e. a, b > 0, a 2 + b 2 = 1. Due to imperfections, the produced state will rather 
be 

Q p , d =p\y)(V\ + (l-p)a, (4) 

where we assume that the noise state a lies inside a separable ball of radius d 
around the totally mixed state, i.e. \\a — 1/4|| < d for some norm. Our aim 
is to provide a local experimental method to tell whether the produced state is 
entangled or not, based on witness operators. 

Since we only want to explain our basic idea, we assume here that we have 
white noise, this means d = 0. The case d > is studied in greater detail 
in [TQ. We would like to stress that our assumption d — is in some sense 
artificial. By this we mean that if d = there is a simple way of determining 
whether q p .q is entangled or not: One can just measure any operator A (which 
fulfills Tr(A) 7^ 4Tr(A\ip) and compute p from the expectation value of 
this operator. With this information the density matrix g Pt Q can be constructed 
and the PPT criterion can be used to decide whether it is entangled or not. 
This is not possible for obvious reasons if d > 0. 

To reach our goal we first have to construct a proper entanglement witness. 
For NPPT entangled states, i.e. entangled states with a non-positive partial 
transpose, the construction of a witness is relatively easy: The partial transpose 
of a projector onto an eigenvector of g Tn with negative eigenvalue does the job. 
Here Tb denotes the partial transposition with respect to subsystem B. 

The state g p $ has one possibly negative eigenvalue, namely A_ = (1— p)/A — 
abp, with the corresponding eigenvector 

|V>_) = -^(|00) - |11)) (5) 

that is independent of a and p. Then Wo = | -0_ )(?/'- | Tb is an entanglement 
witness detecting g Pt o since 

Tr(\^-){^4 TB )Q vfi =TT\^)^-\(Q T p %) = \- <0. (6) 

Having constructed the witness, all that remains is to decompose it. Defining 
|z±) = |0,1), |a; ± ) = (|0)±|l))/ v / 2and 1^) = (|0)±i|l))/ v / 2 we can decompose 
the more general witness \<p)(<f>\ TB for \(j>) = a\QQ) + (3\ll) as follows: 

\4>){cp\ TB = a 2 |00)(00| + /3 2 |ll)(ll| + a /?(|01)<10| + |10)<01|) 

= a 2 \z + z + )(z + z + \ + (3 2 \z- Z -)( Z -z-\ +a(3(\x + x + )(x + x + \ + 
+\x~x~)(x~x~\ - \y + y~)(y + y~\ - \y~y + )(y~y + \) 

= - (l®l + o- z ®a z + {a 2 ~l3 2 ){a z ® 1 + 1® a z ) 

+ 2af3((J x ®<7 X +(Jy® a y )) , (7) 

where the Cj are the Pauli matrices. This way of decomposing a witness can be 
used for higher dimensions and for systems of more than 2 parties by using 

|01)(10| + 1 10) (01 1 = \x + x + )(x + x + \ + \x-x-)(x- X -\- 

-\y + y-)(y + y~\-\y~y + )(y~y + \- (8) 

The local measurement of the general witness of Eq. Q requires three measure- 
ments of Alice and of Bob in the x, y and z direction. This is also true for the 
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special case a = —(3 = l/y/2 corresponding to Wq. It is not possible to evaluate 
the witness with less than three LvNMs: 



Proposition 1. The witness Wo can not be decomposed into less than three 
LvNMs, therefore the decomposition (JJJ) is optimal. 

Proof. The proof was first given in |12| , we repeat it here because we extend it 
later to three qubit systems. Consider a decomposition requiring two measure- 
ments: 



1 

E 



\Al)(Al\®\B}) 



' B j\ 



(9) 



=o 



i,i=0 



where {|^4i)} and {|£?i)} arc orthonormal bases for Ha and Hb, respectively. 
With the help of a Schmidt decomposition as above we can write \tjj)(tp\ TB = 

X« ,=o a i ® a 3 witn 
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Note that the 3x3 sub-matrix in the right bottom corner is of rank 3. Now we 
write any projector on the rhs of Q as a vector in the Bloch sphere: \Aq) (A^ | = 
£$=o s t a i i s represented by the vector s A ° = (1/2, sf , s^, Sg 4 ) and |j4];)(AJ| by 
s A i = (1/2, — sf , — sf, — s A ); \Bq)(Bq\ can be written similarly. If we expand 
the first sum on the rhs of © in the (cr, g) cr, ) basis, the 3x3 sub- matrix in the 
right bottom corner is given by (cq — Cq X — c} + c} 1 )(s^ 1 , s^, S3 ) T (sf , , s^). 
This matrix is of rank one. The corresponding sub-matrix from the second 
LvNM on the rhs of @ is also of rank one and we arrive at a contradiction: No 
matrix of rank three can be written as a sum of two matrices of rank one. □ 
Please note that the decomposition J2J) requires 6 projectors onto product 
vectors (PPV). By applying the method of it is possible to decompose the 
witness using only 5 PPV: 



MM 



^-4^ J2 144X441 - «/3(|oi)(oi| + |io)(io|), (11) 



i=l 



where we have used the definitions 



14) 
14) 
14) 

cos(0) 
sin(0) 



cos(0)|O) +e~ l - sin(0)|l) 
e _< * cos(0)|O) + e^sin(0)|l) 



(12) 



= 14) + 14) 



(13) 



However, with this decomposition the measurement of the witness would require 
4 local correlated measurement settings, hence the two optimization strategies 
are really different. 
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4 Three qubits 



The state space for three qubits has a much richer structure concerning entan- 
glement than the space of two qubits. Let us briefly remind the reader of some 
well known facts about three qubits. We first consider pure states. There are 
two classes of states which are not genuine threepartitc entangled: The fully 
separable states, which can be written as 

\<t>f s )ABC = \a) A ® |/3)b® \l)c, (14) 

and the biseparablc states which can be written as a product state in the bipar- 
tite system, which is created, if two of the three qubits are grouped together to 
one party. One example is 

\<f>bs)A-BC = \ol)a ® \S)BC- (15) 

There are three possibilities of grouping two qubits together, hence there are 
three classes of biseparable states. The genuine threepartitc entangled states 
are the states which are neither fully separable nor biseparable. Given two 
threepartite states, \<f>) and \tfj), one can ask whether it is possible to transform 
\<fi) into \tp) with local operations and classical communication, without requiring 
that this can be done with probability 1 . These operations are called stochastic 
local operations and classical communication (SLOCC). It turns out ^S] that 
\4>) can be transformed into iff there exist operators A, B, C, acting on the 
space of one qubit with 

= A® B® C\<t>). (16) 

Surprisingly, it was proven in |18) that there are two classes of genuine three- 
partite entangled states which cannot be transformed into another by SLOCC. 
One class, the class of GHZ-states can be transformed by SLOCC into 

\GHZ) = l/\/2(|000) + | 111)), (17) 

the other class, the class of W-states can be converted into 

|W) = 1/V3(|100) + |010) + |001). (18) 

Now we can classify the mixed states according to We define a mixed 
state g as fully separable if g can be written as a convex combination of fully 
separable pure states. A state g which is not fully separable is called biseparable 
if it can be written as a convex combination of biseparable pure states. One 
can, of course, define three classes of biseparable mixed states with respect of 
one of the three partitions as well. Finally, g is fully entangled if it is neither 
biseparable nor fully separable. There are again two classes of fully entangled 
mixed states, the W-class and the GHZ-class. g belongs to the W-class, if it can 
be written as a convex combination of pure W- states, and to the GHZ-class 
otherwise. Taking into account that the set of all states is also a convex set, one 
obtains an "onion" -structure. This structure is shown in Fig. 1. 

In the same reference also witnesses for the detection of GHZ-type and W- 
type states have been constructed. Here a want to compute the optimized 
decompositions of these operators. 
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For the GHZ-class a witness operator is given by 
W G hz = ^1-\GHZ)(GHZ\. 



(19) 



If g is a mixed state with Tr(gWoHz) < the state g belongs to the GHZ-class. 
A decomposition of Wghz can be computed with similar methods as for the 
two qubit case, it yields 



Wghz 



5 ■ 1 <g> 1 ® 1 - 1 <g> a z ® a z - u z ® 1 <g> a z - a z ® o z <g> 1 - 
-2 ■ cr^ ® cr^ ® da; + 1/2 • (c^ + dy) ® (a x + a v ) ® (a x + <j v ) + 
+ 1/2 • (a x - o y ) <g> (a x - Oy) <g> (a x - a y 



(20) 



This witness can be measured with four collective measurement settings. Now 
we have to show that this decomposition is optimal. 

Proposition 2. The witness (|19|) can not be measured with three LvNMs, 
i.e. the decomposition H20|) is optimal. 

Proof. The proof is an extension of the two qubit case. First, we write the 
witness in the Oi ® Oj ® ou basis: Wghz = 1/8^ j fc = o AijfcOj ® a j ® °~k, an d 
from (|20|) we obtain: 
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We denote by A^' red the reduced 3x3 matrices that appear when the first 
row and the first column of A^ is dropped: (A^' red ) tJ := {A^) tJ=1 ...^. In 
the same sense one can define a reduced tensor (^f^)i,j,k '■= {^ijk)i,j,k=i,..,3- 
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Let us now investigate what can be achieved with one measurement setting. 
One measurement setting is of the form 

1 

M = Y, c rst \A r )(A r \<E>\B s )(B s \<E>\C t )(C t \ 

r,s,t=0 
3 

= Y, VijkVi ® a 3 ® a k, (21) 

i 7 j : k— 

Defining s A as the Bloch vector of |Ao)(A)| (and similarly s B and s c for \Bo) (Bo\ 
and | Co) (Co |) and using the same argumentation as in the two qubit case, it is 
easy to see that the reduced 3x3x3 tensor fi^ d is given by 

»m= E Crsti-lY+^sfsfsi". (22) 
r,s.t=0 

Therefore for all k the matrices (Myfc)i,j are °f rank one. 

In order to show that Wqhz can not be measured with three measurement 
settings, it suffices to show that it is not possible to find three 3x3 matrices 
Bi,i G {0, ..,2} of rank one such that A^> red , A^' red and A^ 2 ^ red can be rep- 
resented as linear combinations of the Bi . Let us assume the contrary, i. e. that 
we have three Bi. Since the A^' red span a three dimensional subspace in the 
space of all 3 x 3 matrices, the Bi have to be linear independent (as matrices) 
and have to span the same space. That would imply that any of the Bi could be 
written as a linear combination of the A^' red . But a general linear combination 
of the A^' red is of the form: 




This is of rank one if and only if a = /3 = 
the Bi cannot be of rank one and linear 



(3 \ 

a (23) 
7 / 

= 0. Thus, we arrive at a contradiction, 
independent. □ 



For the investigation of W-states two witnesses were constructed in ^Hj • The 
first one is given by 

W W1 = ^1-\W)(W\, (24) 

This witness detects states belonging to the W-class and the GHZ-class, i.e. it's 
expectation value is positive on all biscparable and fully separable states. The 
optimized decomposition is given by 

Wwi = 24 ( 17 " 1 ® 1 ® 1 + 7 " a * ® Uz ® Uz + 

+3- o z ®\®\ + ?>-\®o z ®\ + ?>-\®\®a z + 
+5 • a z <g> a z ® 1 + 5 • a z ® 1 ® a z + 5 • 1 ® a z ® a z + 
-(1 + a z + a x ) ® (1 + a z + a x ) ® (1 + a z + a x ) 
-(1 + a z - a x ) ® (1 + a z - a x ) ® (1 + a z - a x ) 
-(1 + a z + a y ) ® (1 + a z + Oy) ® (1 + o z + Oy) 

-(1 + a z - a y ) ® (1 + a z - a y ) ® (1 + cr, - o-j,)) . (25) 
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Here, only five correlated measurement settings are necessary. This decomposi- 
tion is also optimal: 



Proposition 3. The witness Wwi can not be measured with four measure- 
ment settings, i.e. the decomposition (|25|l is optimal. 

Proof. The strategy of the proof is the same as for the proof of Proposition 2, so 
we can make it short. First one computes Wwi = l/8^j fc = o ^ijkO~i ® <?j ® 
and the corresponding A^' red , i £ 0, .., 3. This time, it turns out that they span 
a four dimensional space. 

Again, it suffices to show that we cannot find four matrices Bi,i £ {0, ..,3} 
of rank one such that A^°^ red , A^' red , A^' red and A^< red can be represented 
as linear combinations of the Bi. Here, the assumption that we have four Bi fails 
due to similar reasons as above: As above, the Bi have to be linear independent 
and it has to be possible to write any of the Bi as a linear combination of the 
A^' red . A general linear combination of the A^'' red is now of the form 
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6 



A = a 7 , (26) 



and this is of rank one if and only ifa = /3 = 7 = 0, hence we arrive at a 
contradiction. □ 



The second witness for W-class states is given by 

W W2 = ^1-\GHZ}{GHZ\. (27) 

This witness can be measured locally with the same decomposition as 120|) 
substracted by 1/4. If —1/4 < Tr(Ww2Q) < 0, g is threepartite entangled, it 
is either a W-state or a GHZ-state. If Tr(Ww2Q) < — 1/4, g is a GHZ-state. It 
also can serve for the detection of states of the type (1 —p)t/8+p\W)(W\, this 
is explained in |19j . 

Let us conclude this section with a remark about the relationship between 
convertability under SLOCC and the number of LvNMs needed for a local mea- 
surement. One may interpret our results for two qubits in the following way: 
A projector |</>)(^>| can be measured with one or three LvNMs, depending on 
whether it is a product state or not. These two classes coincide with the two 
inequivalent (under SLOCC) classes for two qubits (T^|. One may think that 
SLOCC and LvNM are in this way related. Even our results in \F^, which state 
that the number of LvNMs needed strongly depends on the Schmidt rank of 
\4>) for N x iV-systems may support this conjecture, since for bipartite systems 
the Schmidt rank classifies inconvertible sets under SLOCC. But our work on 
three qubits suggests that this coincidence is just by chance. For a general state 
\iPghz) of the GHZ-class there always exists an orthonormal product basis in 
which it can be written as 

H'ghz) = A o |000) + Aie l9 |100) + A 2 |101) + A3 1 110} + A 4 |1H) (28) 

and for a general W-state \tpw) there exists the same description, but with A4 = 
9 = [2n|. If one has an optimized decomposition of a general \4>ghz)(iPghz \ it 
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should be possible to derive a decomposition of |'0iy)('0wH by setting A4 = 9 = 0, 
this decomposition would need less or the same number of LvNMs. In the other 
direction, this means that for a general \ipw) there exists a \ipGHz) which needs 
at least the same number of LvNMs for a local measurement. But as we have 
shown for \W) there also exists a GHZ-state (namely \GHZ)) that needs less 
LvNMs. Hence, the relation between SLOCC and LvNM seems not to be so 
simple, if there is a relation at all. 

5 Conclusion 

We have studied how three qubit entanglement can be investigated with local 
measurements. For this purpose we decomposed already known entanglement 
witnesses into local measurements. We have shown that these decompositions 
are optimal. By this, we have shown that four measurement settings suffice 
for the detection of true threepartite entanglement and especially GHZ-typc 
entanglement. 

We wish to thank Dagmar Brufi, Artur Ekert, Maciej Lewenstein, Chiara 
Macchiavello and Anna Sanpera for helpful discussions. O.G. wishes to thank 
Ofer Biham and Karol Zyczkowski for the discussions in Ustroh, the organizers 
for their work and last but not least for their financial support. This work has 
further been supported by the DFG (Graduicrtcnkollcg 282 and Schwerpunkt 
"Quantcninformationsverarbeitung" ) . 
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